
wo 99/31354 PCT/US98/25866 



Method to Determine Self-Calibrated 



Circumferential Cased Bond Impedance y 



CROSS-REFERENCE TO RELATED APPLICATIONS 
Not Applicable/ 

STATEMENT REGARDING FEDERALLY SPONSORED 
RESEARCH OR DEVELOPMENT 

Not Applicable. 

10 BACKGROUND OF Ttffi INVENTION 

Field of the Invention 

The present invention relates graerally to oil well drilling. More particularly, the present 
invention relates to determining the quality of the bonding between a downhole casing and the 
surrounding formation. Most particularly, the present invention relates to a new system and method 

15 for obtaining an impedance measurement of the cement bonding between a downhole casing and 
the surrounding formation. 
Description of the Related Art 

Although oil well production has been a boon to mankind, with it has come certain 
deleterious side effects. One problem known in oil field drilling is the contamination of fi^sh water 

20 aquifers by hydrocarbon deposits explored by a wellbore. That is, when an oil well borehole is 
drilled through both an aquifer and a hydrocarbon region, there exists a risk that hydrocarbons will 
flow fi-om the hydrocarbon region to the aquifer. This contamination can have a significant 
negative impact on the envhonment and upon drinking water. 

Figure 1 illustrates the occurrence of svich contamination. The surface of the earth 100, a 

25 borehole 110 passing through a rock lithology 160 including an aquifer region 120 and an oil region 
130, a steel casing 140 residing in the borehole, and cement 150 bonding the casing 140 to the 
surrounding formation 160 are pictured. The steel casing 140 is provided around the periphery of 
the borehole 1 10 after drilUng the weU, in part to prevent the collapse of the weU. This steel casing . 
is held in place by special cement 150. If the cement has proper bonding, the risk of oil flow 170 

30 from the hydrocarbon deposits to the aquifer region is reduced significantly. If, however, the casing 
is improperly bonded, the risk that contamination may occur is much greater. 

It has been found that a casing's bonding properties correspond to the impedance of the 
cement bonding the casing to the surrounding formation. The cement impedance indicating proper 
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bonding varies depending upon the type of cement used. Nonetheless, where the cement has an 
impedance of less than 3, there exists poor bonding. At about an impedance of 3, adequate bonding 
has been achieved. Preferable bonding, however, has an impedance on the order of 4 to 5 for soft 
cement, and 5 to 6 for hard cement Therefore, an invention capable of determining the impedance 
of the bonding would be useful to those searching for ways to determine the bonding properties. 
Further, such an invention could be used to evaluate cementing techniques so that more effective 
methods of cementing the casing to the formation may be developed. 

Some efforts have been made to determine by acoustic tools the bonding properties of 
casing. Nonetheless, there still exists a need for an alternative, fast, and reliable device and method 
to determine the quality of the bond between the casing and the formation surrounding the 
borehole. 

SUMNIARY OF THE INVENTION 

A simplified interpretation of the teachings disclosed herein is that the present invention 
finds cement impedance, and therefore the quality of the casing bond, by mathematically fitting 
predicted theoretical reverberation responses with the measured reverberation responses. 
Because the reverberation characteristics of an acoustic wave depend in part on the impedance of 
the cement bonding, the impedance of the cement may be closely approximated. 

Thus, the present invention comprises a combination of features and advantages which 
enable it to overcome various problems of prior devices. The various characteristics described 
above, as well as other features, will be readily apparent to those skilled in the art upon reading the 
following detailed description of the preferred embodiments of the invention, and by referring to 
the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

For a more detailed description of the preferred embodiment of the present invention, 
reference will now be made to the accompanying drawings, wherein: 

Figure 1 illustrates a leak of hydrocarbons to an aquifer; 

Figure 2 illustrates waveform reflection and reverberation; 

Figure 3 is a graph showing a received acoustic waveform; 

Figure 4 is a diagram illustrating the component parts of Figure 3; 

Figure 5 A is a diagram of a subterranean system built in accord with the preferred 
embodiment; 

Figure 5B is a diagram of the above ground system built in accord with the preferred 
embodiment; 

Figure 6A is a flow diagram of the preferred embodiment; 
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Figure 6B is a flow diagram of the preferred embodiment; 
Figure 7A shows a waveform constructed from three digital data points; 
Figure 7B shows the analog waveform corresponding to three data points; 
Figure 8 is a frequency domain illustration of a received waveform; 
5 Figure 9 illustrates the energy deflection of a first gouge in the casing wall; 

Figure 10 illustrates the energy deflection of a second gouge in the casing wall; 
Figure 1 1 illustrates an off-center casing with respect to the wellbore; 
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
Figure 2 illustrates the acoustic wave path and casing reverberations for a downhole 
10 acoustic wave. Shown are acoustic transceiver 200, well fluid 210, steel casing 220, bonding 
cement 230, and a portion of surrounding formation 240. Acoustic transceiver 200 is part of a 
wireUne device (not shown), while well fluid 210 acts as a transmission medium while occupying 
the area between the wireline device and tlie casing, an area referred to as the annulus. Well fluid 
210, casing 220, cement 230, and formation 240 each has its own impedance, labeled Zm, Zs, Zc, 
15 and Zr, respectively. Also shown is acoustic signal 250, including first reflected portion 260, casing 
reverberation portions 270, 271, 272, 273, 274, 275 and cement wave portions 280, 281, 282. 
Reverberations also propagate through the cement due to cement wave portions 280-282, although 

this is not explicitly shown. 

To measure cement impedance, the acoustic transceiver 200 sends out an ultrasonic impulse 

20 250, with a characteristic bandwidth of 195 to 650 kHz, then switches to the receive mode. The 
impulse frequency should be focused on the expected resonance fi«quency of the casing. The 
sound impulse 250 travels through the well fluid 210 and strikes the casing 220. The largest 
portion of the energy of the impulse is reflected back to the transducer as reflected wave 260 while 
a small amount of signal enters the casing as wave 270. When the well fluid 210 is water, the 

25 reflected wave form has an amplitude of about 93% of the initial impulse. The portion of the signal 
that entered the casing is reflected back and forth between the casing/annulus interface and the 
casing/formation interface, as illustrated by wave reverberations 271-275. At each reflection some 
energy is transmitted through the interface, dependent on the acoustic impedance contrast, and is 
either directed back toward the transducer or out into the cement. The signal inside the casing is 

30 quickly dissipated in this manner at a rate directiy dependent on the acoustic impedance of the 
material outside the casing according to the equation: 

Ri = (Zi-Z2)/(Zi+Z2) (1) 
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Where Ri is the reflective coefficient, and Zi and Z2 are the impedances of the materials at the 
interface in question. 

The acoustic transceiver 200, now acting as a receiver or transducer, sees a waveform 
consisting of a loud initial reflection followed by an exponentially decaying signal with peak to 
5 peak times equaling twice the travel time through the casing. Figure 3 illustrates the measured 
acoustic waveform received at the transceiver 200. If time t==0 is the time of generation of the 
acoustic wave at the acoustic transmitter, then the time Ttran represents the transit time (the time 
for the travel of this acoustic wave to the casing and its time to return back to the transceiver). 
As such, the transit time Tiran provides an indication of the downhole distance between the 

10 acoustic transmitter and the casing wall. Also shown in Figure 3 are the Time Offset, Toff, and 
the Resonance Window, Twin, both of whose significance is explained below. 

Figure 4 illustrates the individual waveforms, both first reflection and reverberations, that 
sum to provide the waveform of Figure 3. To find the waveform received by the transducer, each 
reverberation waveform is added to the initial reflection waveform, with each reverberation being 

15 delayed an amount proportional to the width of the casing. Further, because the acoustic transducer 
is not a perfect transmitter, it "rings" somewhat upon the transmission of an acoustic wave. This 
transducer "ringing" also is included in the detected waveforai, and must be accounted for by the 
present invention. 

Because the cement must be properly bonded around the entire periphery of the casing, it is 
20 necessary to determine the impedance of the cement at multiple locations around the circumference 
of the borehole. Therefore, a downhole acoustic transceiver must transmit greater than one shot or 
impulse to ensure that a proper bond exists between the casing and the formation around the entire 
periphery of the borehole. Preferably, an acoustic tool will fire about ICQ shots around the 
circumference of the borehole. Spot size is about .4" wide and 1" high. If water is the drilling 
25 fluid, a stand-off to 1 V-i^ is sufficient, while a more viscous drilling mud requires a smaller stand- 
off. Because each waveform creates its own wave-train, and each wave-tr£un corresponds to the 
cement impedance at a particular point along the casing circumference, a substantial amount of data 
is created for analysis. As such, a mathematically simplified approach is necessary so that cement 
bonding may be obtained real-time. 
30 Figure 5 illustrates a device built in accord with the preferred embodiment. Shown in 

Figure 5 A is acoustic transceiver 200, analog-to-digital converter 500, a tool rotation sensor 504, a 
tool for measuring the speed of sound in mud 508, a processor 510 for recording start time and 
gain, waveform decompression chip 520, and multiplexer 530. Waveform decompression chip 520 
could alternately be part of a processor. Also shown are downhole transmitter 540 connected to 
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multiplexer 530 and telemetry cable 545. Referring now to Figure 5B, at the surface are located 
uphole receiver 550, demultiplexer 560, transmission line 564 carrying tool rotation information to 
processor 590 for a data log 595, transmission line 568 carrying acoustic mud velocity data to 
uphole processor 590, transmission line 570 carrying gain and start time information to uphole 
5 processor 590, and waveform decompression chip 580. Attached to decompression chip 580 is 
processor 590. Processor 590 generates data suitable for a log 595. 

Referring now to both Figs. 5A and 5B, high frequency pulse/echo transducer 200 
circumferentially collects data of casing reflection and reverberation. This acoustic waveform is 
digitized by analog-to-digital converter 500 and sent to processor 510, which detects the first 

10 reflection from the digitized signal. Processor 5 1 0 then computes the relevant start time and transit 
time. Because the total waveform data are greater than the bandwidth capacity of transmission line 
545, digital compression 520 is performed. Suitable compressions include wavelet and ADPCM 
(Adaptive Differential Pulse Code Modulation) techniques, which work well for smoothly 
varying data. To minimize compression error, the waveform data should be fed backward for 

15 \yhen utilizing the ADPCM algorithm. In addition, tool rotation information and tiie speed of 
sound in the transmission fluid are determined. The compressed waveform from digital 
compression chip 520 is then multiplexed 530 witii the tool rotation information and the 
measured speed of sound in the transmission medium (e.g. mud). Downhole transmitter 540 
sends this multiplexed data to the surface. Sending the data to the surface allows processing by 

20 faster, more sophisticated machinery. 

This mvdtiplexed data is received by uphole receiver 550 and is separated into component 
parts by demultiplexer 560. Waveform decompression chip 580 provides the reconstructed 
waveform to processor 590, which also receives start time and mud speed information. Upon the 
determination of the cement impedance, processor 590 combines this with tool rotation 

25 information and creates a log 595. 

Figure 6 illustrates a preferred method for the present invention. At step 600, ah observed 
waveform is provided downhole for processing. The waveform's transit time (Tiran) is obtained at 
step 610, as well as the time windows Toff and Twin. The definition of transit time was explained 
above with reference to Figure 3 and may be easily measured by a fust reflection detector portion 

30 of processor 510. Totr and Twin are then selected to obtain a time window Twin that contains 
reliable reverberation information. Tom measured from the time of receipt for the initial 
reflection, is a time window that encompasses the initial reflection. As such, its duration is 
dependent upon the frequency of tiie acoustic impulse transmitted by acoustic transceiver 200 
and the nature of the drilling fluid. Toff also must account for error introduced because of the 
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real-world shortcomings of the acoustic transducer (transducer "ringing")^ and thus Toff must be 
slightly longer than if chosen theoretically. Nonetheless, Toff is about 15 microseconds. Twin is 
juxtaposed with Toff and is a time window of interest because Twin contains reverberation 
information uncontaminated by the first reflection. The duration of Twin should be brief enough 
5 so that noise and reverberations occurring in the cement 230 do not make unreliable the received 
casing reverberation waveforms. Nonetheless, so that a reliable wave train containing sufficient 
data is obtained. Twin preferably includes at least four reverberations. Thus, Twin is about 12-8 
microseconds. 

At step 620, the peak amplitude of the initial reflection (contained within Toff)' is 
10 determined. One complicating factor when determining the value of the peak amplitude is the 
digitized nature of the waveform signal received uphole. Further, the waveform data is digitized 
coarsely (i.e. the digital data points are spaced out widely). As such, undesired error is 
introduced. To correct the error quickly and accurately with a minimum of processing, the 
present invention employs a quadratic interpolation approach to approximate closely the 
15 amplitude peak of the initial reflection. 

Figure 7 A shows an absolute value digital waveform 730 modeled on three digitad data 
points a, b, and c, labeled 700, 710, and 720, respectively. Quadratic curve 740 including a 
maximum amplitude point labeled X 745, are also shown. Figure 7B shows the absolute value 
analog waveform 750 as detected downhole (corresponding to digital waveform 730). Also 
20 illustrated are points a, b, and c, in addition to point p 760, the peak amplitude point of the analog 
waveform. 

Referring back to Figure 7 A, points a 700, b 710, and c 720, are transmitted uphole as 
part of the digitized and compressed waveform, and these points are then used to construct a 
waveform 730. However, waveform 730, based on digital data points, is not an exact replica of 
25 the analog waveform 740. As shown in Figure 7B, the absolute maximum amplitude point p of 
the analog waveform 740 is greater than the amplitude corresponding to point b 710. This effect 
becomes more pronounced as the digital data sampling becomes more coarse. To account for 
this error, the present invention uses a quadratic approach based on points a, b, and c to find the 
actual analog peak. Referring now to both Figures 7A and 7B, the peak 745 of the quadratic curve 
30 740 closely approximates the actual peak of the waveform 750, point p 760. Thus, a quadratic 
approximation based on points a, b, and c may be used uphole to approximate the actual peak of 
the anedog waveform 750. 

To approximate the actual peak amplitude with a great degree of precision, three adjacent 
absolute peaks (a, b, c) must be found, where the magnitude of a is less than b, which in turn is 
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greater then c (a<b>c). Then, points a, b, and c are fit to a quadratic curve 740 having the 
relationship: 

Y = A + Bx + Cx^ (2) 

5 

By fitting this equation to the known points a, b and c as shown in Figure 7, the coefficients A, 
B, and C may be found. The amplitude (y-axis) and time (x-axis) of peak 745 may then be found 
from the quadratic curve as follows: 

10 Amplitude of peak =A-B^/4C 

Time of peak = -B/2C 

' Referring again to Figure 6 A, the casing thickness is found at step 630. Casing thickness 
may be found according to the relation: 

15 

CT = Vs/(2fR) (3) 

Where, 

Ct = casing thickness 
20 Vs = the velocity of sound in steel 

fR = resonance frequency of the casing 

The velocity of sound in steel, Vs, is a known quantity. The resonance frequency of the 
casing may be obtained as follows. First, the reverberation waveform defined by Twin is 
transformed to the frequency domain by use of DFT (Discrete Fourier Transform), Figure 8 

25 illustrates a frequency domain graph of such a reverberation waveform. As is well known, the 
DFT embodies the proposition that a waveform may be represented by an infinite number of 
sinusoids at varying frequencies. Because resonance frequency is defined as the frequency at 
which maximum energy is reflected back by the casing, the resonance frequency corresponds to 
frequency at which peak amplitude is obtained. This frequency is labeled fa. Frequency fR 

30 corresponds to an amplitude point labeled XfR. Also shown are frequencies fa, fb, and fc To fmd 
the peak amplitude point Xm, an initial guess of resonance frequency is used to establish a first 
amplitude, Xfa. By incrementing (if Xfa < Xfb) or decrementing (if Xfa > Xn,) the guess frequency to 
fb, a second amplitude value Xfb is obtained. Upon obtaining a third amplitude Xfc satisfying the 
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relationship Xfa<Xfb>Xfc, a quadratic approximation may once again be employed to find the 
amplitude peak XfR and corresponding resonance frequency fR. 

The accurate determination of resonance frequency and peak amplitude, combined with the 
measured digital data, allows an accurate reconstruction uphole of the analog reflection signal 260 
contained in Toff. The accurate reconstruction of the reverberation waveform contained in Twin is 
not as problematic. 

Referring back to Figure 6A, proper modeling applied to the first reflection signal 260, as 
defined by Totr, gives a theoretical prediction of what the reverberation waveform contained in Twin 
should look like. To accomplish this, at step 640 the first reflection signal is transformed by Fast 
Fourier Transform (FFT) into its frequency domain equivalent. This yields S(co). Because the 
modeling is done in the frequency domain, amplitude and phase errors are eliminated. This error 
elimination simplifies mathematical processing (and hence faster processing is obtained). 

Altemately, instead of transforming each first reflection individually, to simplify and speed 
mathematical processing, the first reflections may first be averaged and the result transformed at 
step 640 by FFT processing into the frequency domain to yield S(co). Further, averaging of the first 
reflection is preferred because, due to the same transducer collecting all the data along the 
circumference and the majority of the energy being in the first reflection, most of the transducer 
properties are in the initial reflection. Thus, the average of the first reflections along the 
circumference will give a better estimate of the transducer properties, and thus a better input to the 
theoretical modeling. A most reliable first reflection average may be obtained by discarding first 
reflections that have amplitudes above or below a preset deviation from a moving average of 
preceding first reflections (i.e. by discarding those first reflection waveforms that correspond to 
some type of defect in the casing). 

At step 650, a theoretical prediction of the reverberation waves is obtained by convolution 
of the frequency-domain first reflection signal, S(co), with a fi^quency-domain theoretical response 
equation, R(co). Assuming a flat casing, the theoretical frequency domain response may be 
modeled by the following: 



4Z,Z, Z,-Z, 



e 




(4) 



Z„+Z, Z,+Z, 



1- 



Where 
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R((0) = the reflection coefiFicient for angular frequency to 

Zm, Zs, Zc = impedances for mud, steel casing and cement behind casing 

Vs = the speed ofsound in the steel casuig 

5 Although the casing has a degree of curvature (because it is cylindrical), it may accurately 

be presumed as a flat surface locally in view of the comparatively short wavelengths of the 
acoustic signal. This is because the frequency and shape of the transducer, combined vidth the 
offset between the transducer and the casing, ensure that the acoustic beam does not spread 
significantly. For smaller casings, more severe casing curvature may be compensated for by a 

10 smaller offset between the transducer and the casing wall, thereby precluding appreciable beam 
spreading. Only where there are streaks, roughness, or gouges in the casing must the present 
invention normalize the received waveform, as described below with respect to step 680. 

For equation (4), Vs is a well known quantity and Ct has been obtained for a particular 
waveform at step 630. The use of an average of casing thickness, Cxav, may instead be used at 

15 step 650. The use of an average of casing thickness helps simplify and speed mathematical 
processing, making obtainable real time derivations. The use of Cxav is also preferred because 
casing thickness does not vary significanUy along an azimuth. As such, average casing thickness 
is the best estimate available to calculate the coefficients ao, bo, co, and do, explained below. Zm 
may be obtained upon determining Vm, the speed of sound in mud, by use of the following 

20 relationship: 

Zm = Vmp (5) 

Where, 

25 Zm = impedance ofthe transmission medium (mud) 

Vm = velocity of an acoustic wave in the transmission medium 
p = density ofthe transmission medium 
Methods to determine p by measurement, and Vm by tool measurement 508, are well known. 

At step 660, these results are transformed back into the time domain by use of an Diverse 
30 Fast Fourier Transform (IFFT). As such, step 660 provides the theoretical reverberation response 
for the observed initial reflection waveform(s) in the time domain, the tiieoretical resonance sum 
being represented with the upper case notation ■ ™^ theoretical reverberation response, 

obtained at step 660, is used as a yardstick for the observed reverberation waveform contained in 

Twin. 
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Next, four unknowns are obtained at step 665. As an initial matter, for a reverberation 
window of interest. Twin, the natural log of the sum of the reverberation waveform amplitude varies 
linearly with cement impedance. That is, a linear relationship between cement impedance and Sw 
may be expressed as: 

5 

Zc = A + B In (Sw) (6) 
where Sw is the sum of the reverberation waveform amplitudes and has the form: 

10 Sw= XKI 

the lower case Xi being the amplitude at any given point in the observed reverberation waveform 
contained in Twin- If, 

15 A = ao + boCt (8) 

B = co + doCt (9) 

then equation (6) may be expressed as: 

20 Zc = ao + bo Ct + co In(Sw) + do Ci In(Sw) (10) 

Where Zc is the cement impedance behind the casing. 

Substituting the theoretical resonance sum in lieu of Sw, the right side of equation 

(10) contains Ci and four unknowns: ao, bo, co, and do. At step 665, these four coefficients are 

25 determined by proper substitution and use of the well known least square approximation 
technique. More specifically, the determination of four unknowns, ao, bo, co, and do, requires at 
least four equations. A higher order approximation for the values of ap, bo, co, and do may be 
obtained by the use of more than four equations. For example, the values 1, 4, and 8 may 
variously be introduced as values for Zc, and Cxav +A dCrav as established at step 630 may be 

30 substituted for three values of Ci. The magnitude of dCjav is about 0.025 inches. These 
substitutions, then, yield nine equations. Use of the well known least square approximation 
techniques yields accurate estimates of ao, bo, co, and do. 

Next, the actual resonance sum Sw is obtained at step 670. As defined above: 
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Sw= SKI 

where the lower case Xi is the observed/measured amplitude at any given point in the 
reverberation waveform. 

However, scratches, roughness, and gouges can affect the detected reverberation 
5 waveform. For example, a gouge such as that shown in Figure 9 can reduce the amount of 
energy detected by the transducer 200. A gouge such as that shown in Figure 10 can increase the 
amount of energy detected by the transducer 200. Such errors could lead to an indication of a 
free pipe when, in truth, a good bond is present, or vice versa. 

Other conditions also lead to similarly misleading initial reflections, such as a non- 
10 circular casing (not shown), or the tool being off-center with respect to the borehole, shown in 
Figure 11. The circular geometry or lack thereof for the casing may be established by a 
comparison of transit times for the initial reflection, Tiran. 

Figure 11 includes a weilbore 1100, casing 1120 in wellbore 1100, a well bore center 
point 1 105, and a casing center point 1125. Also shown are normal lines 1130, 1 140 and 1 150; 
15 normal line 1 130 originates from wellbore center point 1 105, normal line 1140 originates from 
casing center point 1 125, and normal line 1 150 travels through both center points 1 105 and 1 125. 
Graph 1160 includes amplitude and azimuth axis. Lines 1170 and 1175 correlate the two 
normal coincidence angles and the two maximum amplitude points Amax- 

Where the casing is off-center with respect to the borehole, the initial reflection amplitude 
20 will be highest when both normal lines coincide, such as at normal line 1150. Otherwise, the 
amplimde is decreased. As such, a sine wave amplitude pattern emerges, with the amplitude 
peaks occurring where the normal lines coincide. This phenomenon may be used to identify the 
off-center casing with respect to the wellbore. Thus, the present invention corrects for the 
presence of irregularities in and around the casing to obtain optimum impedance measurement 
25 accuracy. 

The resonance sum is corrected to account for defects at step 680. This is complicated 
because based upon a single observation, it is extremely difficult to determine whether the casing 
has a defect. Further, the presence of a casing defect is most apparent from examination of the 
initial reflection. As such, to help determine whether a defect exists, the initial reflection 
30 waveform is compared to a moving average of a number of previous reflection amplitude 
measurements. In other words, a moving average filter using a limited number (e.g. five) of 
immediately previous measurements is compared to the received initial reflection. If the received 
waveform amplimde differs more than a preset amount from this moving average, a correction 
for irregular peak amplitude is necessary. The initial reflection waveform is corrected by 
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normalizing its amplitude to that of the moving average. This yields the magnitude of the 
correction for the initial reflection waveform. The amount of correction necessary for the 
resonance sum may be established based on a linear relation. That is, the amplitude of a 
resonance sum is directly proportional to the peak amplitude of the first reflection variation with 

5 respect to the average peak amplitude. As can be readily understood, as the initial waveform 
amplitude is normalized in a particular direction, the resonance sum is likewise adjusted. For 
example, if the first reflection peak amplitude is 'A the average peak amplitude, the resonance 
sum must be corrected by doubling its original value. So that subsequent resonance sums are 
accurately corrected, any significantly different amplitude initial reflection waveform should be 

10 left out of the computations for subsequent moving averages. Lastly, at step 690, Sw is 
substituted into equation (10), and cement impedance Zc is determined for any point along the 
circumference of the casing. 

While preferred embodiments of this invention have been shown and described, 
modifications thereof can be made by one skilled in the art without departing fi-om the spirit or 

15 teaching of this invention. The embodiments described herein are exemplary only and are not 
limiting. Many variations and modifications of the system and apparatus are possible and are 
within the scope of the invention. Accordingly, the scope of protection is not Umited to the 
embodiments described herein, but is only limited by the claims that follow, the scope of which 
shall include all equivalents of the subject matter of the claims. 

20 
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CLAIMS 

WHAT IS CLAIMED IS: 

1 . A method for determining casing bond impedance, comprising: 

(a) generating an acoustic signal within a borehole casing, said casing having a 
bond impedance; 

(b) receiving reflections of said acoustic signal from said casing, said reflections 
comprising an initial reflection portion and a reverberation portion; and 

(c) analyzing said reverberation portion to determine said bond impedance. 
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